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A  curvilinear  regression  model  is  treated  by 
linear  programming  methods,  so  as  to  require 
only  mild  restrictions  on  the  nature  of  the 
non-linearities. 


1.  Introduction 

Consider  the  regression  model 


(1) 


+  error 


where  y  is  the  dependent  variable,  and  Xj  ,  J  =  1,2,. ..,p  ,  are  the 
independent  variables.  Assiirae  there  are  available  k  sets  of  observa¬ 
tions  on  the  variables,  viz.,  y^^  and  ^  i  =  l,2,...,k  ,  and  that 

we  wish  to  estimate  the  functions  f  (x  )  from  these  data.  We  do  not 

U  J 


One  of  the 
i  =  1,2,. ..,k 


X  may  represent  a  dummy  variable. 


Xij  ,  -  1 


impose  the  restriction  that  every  fj(Xj)  must  be  linear,  and  thus  our 
model  may  be  categorized  as  curvilinear  regression.  The  usual  statis¬ 
tical  technique  for  such  models  is  to  specify  a  mathematical  form  for 

2 

each  of  the  functions  alogx^  ,  and 


% 


and 


then  to  estimate  the  parameters  In  the  specified  forms  according  to  a 

least  squares  criterion.  An  alternative  approach,  given  In  [4],  employs 

criteria  that  permit  the  use  of  linear  programming  models  for  obtaining 

the  estimated  parameter  values.  In  this  paper  we  do  not  go  so  far  as 

specifying  a  mathematical  form  for  the  f  (x.)  ;  rather  we  Impose  only 

J  J 

mild  restrictions  on  the  behavior  of  the  functions;  for  example,  a  given 

f.(x  )  must  be  monotonlcally  Increasing.  Utilizing  a  linear  programming 
J  J 

approach,  we  are  able  to  estimate  the  functions  f  (x. )  given  only  the 

J  J 

minimal  assumptions  as  to  the  required  shape.  In  Section  2  we  adopt 
"minimal  sum  of  absolute  deviations"  as  the  criterion  of  best  fit;  In 
Section  5  ''s  review  alternative  criteria. 

2 .  Basic  Formulation 

To  begin,  suppose  that  for  each  J  the  values  of  x^^j  are  distinct. 


values  for  the  f . , 
^  J 


CXir  procedure  is  to  find 


For  each  observation  y^  ,  we  have  an  Error  Constraint, 


(2) 


^  ^ij  ^11 
J=1 


"21  “  ^1  >  ^  -  l;2,...,k  , 


"li  ^  O  '  "21  ^  0  • 

As  Is  explained  in  detail  in  [4],  and  are  interpreted  as 

vertical  deviations  "above"  and  "below"  the  fitted  regression  function 
for  the  1-th  set  of  observations;  i.e.,  the  sum  +  eg^^^  represents 


the  absolute  deviation  between  the  estimated  fit  for 
and  y^  .  Our  criterion  of  best  fit  is  a  minimal  sum  of 
the  absolute  deviations. 


J=1 
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Next  we  impose  Weak  Curvature  Constraints  on  the  .  To 

illustrate,  we  suppose  for  J  »  1,  <  x^j^  <  . . .  <  Xj^  ,  and  we 

postulate  that  f^(Xj^)  is  a  monotonic  increasing  function.  Bien  the 
Weak  Curvature  Constraints  for  the  set  of  f^^^^  are 

(3a)  ^il  5  ^i+1,1  "  i  =  1,2,.  ..,k-l  . 


In  other  words,  we  Impose  by  (3a)  a  complete  ordering  on  the  f  .  By 

i  J 

reversing  the  Inequality  in  (3a),  we  thereby  would  impose  the  restriction 
that  f^(x^)  be  monotonically  decreasing.  If  instead  we  hypothesize 
that  f^(x^)  is  unlmodal,  the  peak  occurring,  for  example,  at  a  value  of 
X  between  x^^  and  X|^j^  ,  the  Weak  Curvatvire  Constraints  for  the  set 
of  fj^j^  are 


(3b)  i  -  1>2;  ^  i  =  U,5,  ,..,k-l. 

Notice  that  (3b)  depends  on  being  able  to  specify  that  the  modal  vsCLue 
occurs  for  an  x  between  two  given  values  of  x^^^^  .  The  estimation  of 
does  not  yield  the  modal  value  unless  x  equals  one  of  the  given 

Xii  . 

Momentary  reflection  is  sufficient  to  demonstrate  that  the  fore¬ 
going  analysis  is  modified  in  a  straightforward  manner  if,  for  a  given 
j  ,  the  x^j  values  are  not  ordered  monotonically.  For  example,  if 
ordering  the  x^^  yields  x^^  <  x^^  <  Xg^  <  ,  and  we 

postulate  that  ^^^(^q^)  is  a  monotonically  increasing  function,  then 
instead  of  the  Weak  Curvature  Constraints  (3a),  we  have 

^31  =  ^11  =  ^81  =  ^91  =  ^21  ■ ’ ■  ■ 


-3- 


similarly.  It  Is  possible  by  a  set  of  Inequality  constraints  to 

Impose  partial  orderings  on  the  f .  .  to  produce  blmodal  functions, 

^  J 

fiinctlons  having  both  a  single  relative  minimum  and  a  single  relative 
maxlmvun,  etc . 

Up  to  this  point,  we  have  assumed  that  for  each  J  the  vaG.ues  of 

the  X.  are  distinct.  If,  for  a  particular  J  ,  two  or  more  values 
^  J 

of  X  are  Identical,  then  the  corresponding  f  must  be  equal, 
ij  IJ 

Consequently  for  each  J  It  Is  only  necessary  to  define  a  variable 

fj^j  In  (2)  wherever  Identical  values  for  occur.  Likewise  In  (3) 

It  Is  only  necessary  to  define  Inequalities  for  the  distinct  f .  . 

^  J 

Cie  basic  linear  programming  model  may  now  be  stated  as  finding  values 

for  the  f  which 
^  J 

k  k 

(U)  minimize  c  +  JZ 

1.1  1=1 

subject  to  the  linear  constraints  (2)  and  (3). 

We  complete  our  discussion  of  the  basic  formulation  by  determining 

the  nizmber  of  restrictive  relations  and  variables  that  comprise  the  model. 

For  the  sake  of  definiteness,  suppose  that  all  the  functions  f.(x  ) 

<J  J 

are  postulated  to  be  monotonic.  Let  N  be  the  total  number  of  distinct 
f^j  to  be  Included  in  the  constraints.  Hien  (2)  and  (3)  Imply 


Specifically,  for  each  j  ,  let  n  be  the  total  number  of  distinct 

0 

values  of  ;  then  N  =  n 

ij  J=1  ^ 


Number  of  Equalities  =  k  (Error  Constraints) 

Number  of  Inequalities  =  N-p  (Weak  Curvature  Constraints) 
Number  of  Variables  =  N  +  2k  ^21 
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If  slack  variables  are  added  to  the  Weak  Curvature  Constraints  to  convert 


them  to  equalities,  we  have 


Number  of  Equalities  =  k  +  N  -  p 
Number  of  Variables  *  2(N  +  k)  -  p  . 

Although  In  some  real  regression  applications  these  numbers  need  not  be 
negligible,  the  method  remains  feasible  since  current  linear  programming 
computer  codes  are  capable  of  solving  several  hundred  equations  In  many 
hundreds  of  unknowns. 

In  our  discussion  so  far  we  have  restricted  only  and  to 

be  non-negative.  In  the  next  section  we  discuss  placing  analogous 
restrictions  on  the  f^^j  ,  but  here  we  mention  that  If  an  f^^j  Is  un¬ 
restricted  in  sign,  then  a  common  way  of  treating  it  within  the  usual 
format  of  a  linear  programming  model  having  all  non-negative  variables 
is  to  let  f.  .  =  f^.  -  f7,  ,  f^,  >  0  ,  fj,  >  0  ,  i.e.,  to  substitute - 

into  the  model  restrictions  in  place  of  f^j  ,  the  difference  between 

two  non-negative  variables.  Ihis  device  increases  the  number  of  variables 
but  has  no  effect  on  the  number  of  constraints.  We  return  to  the  question 
of  model  size  in  Section  U,  where  we  discuss  the  dual  formulation. 

3-  Additional  Constraints 

We  may,  if  we  wish.  Impose  stronger  properties  on  the  estimated 
values  for  the  f  (x  )  functions.  We  illustrate  several  such  possibilities. 
We  may  add  Positioning  Constraints  such  as,  for  each  j  ,  some  or  all  of 
the  f^j  must  be  non-negative,  the  smallest  f^^j  must  be  no  smaller 

than  Lj  ,  and  the  largest  f^^j  must  be  no  larger  than  Uj  .  These 

restrictions  are  all  representable  as  linear  Inequalities. 
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We  may  apply  Strong  Curvature  Constraints  Instead  of,  or  In  addition 


i  t 

i 


(S)-  For  example.  If  the  estimate  for  strictly 

monotonlcally  Increasing  and  the  minimal  allowable  difference  between  two 
adjacent  values  of  f^^^  Is  to  be  6^^^  ,  then  (3a)  would  be  modified  to  the 


linear  relation 


(3<!) 


“■ii  *■  '‘ii  S  '■1*1,1 


similarly,  if  there  is  to  be  a  maximal  allowable  difference  between 

two  subsequent  values  of  f^^^  ,  then  we  would  add  the  linear  constraints 

^3*^^  ^11  *  ^1+1,1  ■ 

The  relations  (3c)  siid.  (3d)  imply  absolute  differences  between  adjacent 
values  for  f^^^  .  It  Is  also  possible  to  Impose  constraints  on  the 
relative  differences  between  values  for  ;  to  illustrate,  we  may 

require 


<*■21  -  ^l)  S  ■  ^2l) 


(3*) 


(^21  -  '■11)  S  ^^Sl  -  f21> 

where  a  and  p  are  given  positive  constants. 

Finally,  it  is  possible  to  specify  an  exact  mathematical  form  for 

any  particular  f  (x  )  .  For  example,  if  we  postulate  that 

J  o 
o 

f^(x^)  =  ax^  ,  then  we  modify  (2)  to 
(5) 


^li  ■  ^21  -  ^  -  1^2,..., k  , 


and  eliminate  the  Curvature  Constraints  for  x^^^  ;  the  linear  programming 
model  treats  a  as  a  variable  in  the  optimization  process. 


( 
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I*-.  Dual  Formulation 


As  we  have  demonstrated  earlier  [4],  there  may  he  computational 
efficiencies  arising  from  solving  the  dual  linear  programming  problem  [3] 
in  the  process  of  optimizing  the  dual  problem,  an  optimum  solution  to  the 
original  regression  model  is  automatically  generated  [14-].  We  demonstrate 
the  argument  with  reference  to  the  example  at  the  end  of  Section  2, 
postulating  further  the  restrictions  f^j  >  0  .  Let  w^  denote  the  dual 
variables, where,  in  the  primal  model,  i  =  1,2, ...,k.  correspond  to  the 
Error  Constraints  and  i  =  k+l,k+2, . . . ,k+N-p  correspond  to  the  Weak 
Curvature  Constraints.  Then  the  dual  relations  associated  with  the  error 
variables  and  are 

(6a)  w^  <  1  ,  -w^  <1  ,  1=1, 2,. ..,k 

We  make  a  transformation  of  variables 

(7)  =  w^  +  1  ,  i=l,2,...,k  , 

and  substitute  (z^^  -  l)  for  w^  first  into  the  N  dual  relations 
associated  with  the  variable  f^^  ,  and  then  into  (6a),  the  latter 
yielding  a  set  of  upper  bounds 

(6b)  °  5  ^i  ^  2  ,  i  =  1,2,.  .  .,k  . 

As  a  result,  we  have  in  the  dual  model 

Number  of  principal  inequalities  =  N  (corresponding  to  the  primal 

f^j  variables) 

Number  of  upper  bound  inequalities  =  k 
Number  of  variables  =  k  +  N  -  p 

If  slack  variables  are  added  to  the  principal  inequalities,  we  have 
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Number  of  principal  equalities  =  N 
Number  of  upper  bound  Inequalities  =  k 
Number  of  variables  *  2N  +  k  -  p 

Therefore  if  the  dual  problem  is  solved  directly  instead  of  the  primal 
model,  a  reduction  of  (k-p)  equations  and  k  variables  occurs  at  the 
extra  computational  expense  of  handling  k  bounded  variables.  There 
are  efficient  algorithms  for  solving  bounded  variables  problems  (see 
[2]  and  [3])  and  a  number  of  computer  codes  are  available.  If  additional 
constraints  are  appended  in  the  primal  problem,  such  as  those  suggested 
in  Section  3,  then  utilizing  the  dual  problem  becomes  a  more  attractive 
approach  since  these  relations  increase  only  the  number  of  dual 
variables  and  not  the  nxmiber  of  dual  restrictions. 


5.  Alternative  Objective  Functions 

Instead  of  employing  a  minimal  sum  of  absolute  deviations  as  a 
criterion  of  best  fit,  it  is  possible  to  use  a  Chebyshev  or  a  least 
squares  criterion.  Under  a  Chebychev  criterion,  we  seek  to  find 
values  for  the  f^j  which  minimize  the  maximum  absolute  error  between 
the  fitted  relation  and  y  .  To  accomplish  this  aim,  we  have  in  place 
of  (2)  pairs  of  Error  Constraints 


-  foj  *  '  i  -yi 


J=1 


(8) 


'  i  Ji 


i  —  1,2, . . . ,k  , 


J=1 


€  >  O 


and  in  place  of  (*^),  the  objective  function 
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(9)  minimize  c 

Con'binuing  with  our  example  wherein  there  are  N  distinct  , 

and  f . ,  >  0  ,  we  have  in  the  primal  model 

i  J  ® 

Number  of  inequalities  .  2k  +  N  -  p  (Error  and  Weak  Curvature 

Constraints) 

Number  of  variables  =  N  +  1  • 

In  the  formulation  the  above  values  for  the  number  of  inequalities 

and  variables  are  interchanged,  leading  to  the  conclusion  that  the  dual 
is  in  general  computationally  more  efficient  to  solve. 

Under  a  least  squares  criterion,  the  only  change  in  the  basic 
formulation  of  Section  2  is  the  objective  function,  which  is  then 
stated  as  follows; 


2  ^  2 

(10)  minimize  (  YZ.  ^2i  ^  ’ 

i=l  i*l 

This  problem  is  a  quadratic  programming  model,  and  can  be  solved  by  a 
method  such  as  that  in  [5],  for  which  computer  routines  are  available. 
It  is  also  possible  to  employ  a  linear  programming  version  if  the 
squared  functions  in  (lO)  are  represented  in  terms  of  piecewise  linear 
approximations  ( 1  ] . 


6.  Concluding  Remarks 

The  model  in  this  paper  permits  a  large  degree  of  flexibility  in 

* 

estimating  non-linear  regression  fiuactions.  If  the  underlying  purpose 


We  point  out  that  the  approach  is  well  suited  to  situations  in  which 


the 


'ij 


represent  rank  order  measurements. 


for  finding  the  functions  is  entirely  pragmatic,  then  one  may  utilize 
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the  solution  values  for  the  f, 


without  further  consideration.  To  the 


1 


ij 

extent  that  the  ultimate  purpose  may  be  to  discover  simple  "regression 


laws,”  then,  once  the  values  for  the  f^j  are  found  and  plotted  against 

X  ,  one  may  wish  to  postulate  a  particular  mathematical  form  for  each 
^  J 

f.(x  )  function,  and  estimate  the  parameters  in  the  mathematical  forms 
J  J 

either  by  least  squares  or  by  linear  programming  methods. 
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Chicago  1,  lUinela  1 

Stvarlntandaet 

0.  B.  laval  Peatgradoata  School 
Attn:  Xdbrary 

Monterey,  Callfemia  1 

Sleetroaie  Convutar  Dlvlaieo 

Coda  260 

BuraaM  of  Shlpa 

Oapartaant  of  tha  Bavy 

WhahlagtoB  25,  D>  C,  2 

the  RAID  Corporation 
1700  IMln  Straat 

Santa  Ncalea,  California  1 

Miiasv— Controla,  laf.  Library 
Baalagtcn  Rand, 

Dlv.  of  Sparry  Rand  Corp. 

313  Pourtb  Avataia 

Bar  Yof^  10,  1.  Y.  1 

Unlvaraity  of  Warhlagtrin 

Inat.  for  BeonoMle  Raaaareh 

Seattle  3,  khahU^toa  1 

Unlvaraity  of  California 
Bureau  of  Bualneaa  and 
BconcBile  Raaaareh 

Berkeley  k,  Callfcnla  1 

Ifr.  Ads  Arad 

Applied  Statlatlca  Qro«A 

Dapartaant  of  Mmagsant 

Stavaaa  Xaatltvta  of  IhchBolegy 

Bobokan,  law  Jaraay  1 

Dr.  8.  O.  Allan 

Stanford  Raaaareb  Xnatltiita 

Manlo  Park,  California  1 

Prefasa%tf  laland  Baeb 
Dapartaant  of  Bcoooalea 
Camagla  Znatltvta  of  Tachnolcv 
PlttWbvgh  13,  PanBaylvanla  1 

Profeaeor  B.  P.  Rarltanliarti 
Dapartaast  of  Mathaaatlea 
tkilvaral^  of  California 

Loa  Angalae  2%,  California  1 

Dr.  Martla  J.  Acknan 
>partai  nt  of  iceasics 
M-ew«  Qhlvaralty 

Provldaoea  12,  Mwada  Xalaad  1 

Dr.  Richard  RallaaB 
Iha  im  Cue  pui  ntlen 
1700  Mala  Street 

Santa  Monica,  Oallfomla  1 


Profaaaor  Ms  1.  Bos 

Collapa  of  ■aglnaaa  AiMlalatmtloa 

flpiaauga  Dhivnralty 

In' ar  nan  10,  law  Yert  1 

fiaan  L.  N.  K.  Bsalter 

School  of  MMlBoarlM 

Dhlvaratty  of  California 

Leo  AngalM  2k,  California  1 

Profaaaer  Jsna  Bole# 

Dtalvaralty  of  California 

iplenltaral  Mvarlaant  Statls 

Bnrkalay  k,  California  1 

Profeaeor  0.  I.  Mrownlaa 
Dapartaant  ef  Beans! ea 
Unlvaraity  ef  Mlnaaaota 
MisaapoUa  Ik,  Nlnnaaota  1 

Prefaeaor  8.  8.  Calma 
■aad.  Dept,  ef  IMthaamtlea 
Unlvaraity  of  XUlnola 

Urbana,  llllaela  1 

Mr.  Bssl  L.  Carrell 
7022  Yokalen  Drive 

Dalian,  Vaxaa  1 

Profeaeor  A.  Chamaa 
The  Vachnelegieal  Inatltuta 
Borthwastarn  Unlvaraity 

Bvmnator,  XUlaela  1 

Profaaaer  John  B.  Chlpnan 
Ihpartnant  of  leeasiea 
unlvaraity  ef  Mlnnaaeta 

Minnaavolla  Ik,  Mlnmaota  1 

Profeaeor  Carl  Crlat 
Bahavioral  Selaneaa 
202  Jslpare  Barra 

Stanford,  California  1 

Prefaaaer  Randolph  Chnrcb 

U>  8.  lavnJ  Peatgraduata  Bebool 

Montaray,  California  1 

Profeaeor  C.  V.  Cbarefcnan 

School  ef  Bnainaae  Adnlnlatratls 

UnivaraltT  ef  Callfea'nla 

Barkalay  k,  California  I 

Profeaeor  V.  V.  Cooper 

Director  of  Raaaareh 

Dapartaant  ef  Eeooealca 

Camagla  Inatltuta  of  Taebnalogy 

PittA«a*ih  13,  Pannaylvanta  1 

Profasaor  Oarald  Dehrau 

Covlea  Potsidatioo  for  Raaaareh  In  lea. 

Yale  Station,  Boa  2123 

law  Bavan,  Connaetleut  1 

Profaaaer  Robert  Dorfhaa 
Dapartaant  of  Bcenealea 
Harvard  Unlvaraity 

Caabrldga  38,  Maaanchunatta  1 

1^.  Loola  Xkron 
IMnaganant.  fnfftnaarlTif  Staff 
Baras  of  Tarda  aad  Doeka 
lavy  Dapartast 

Mhahlngton  25,  D.  C.  1 

Raar  Mn.  B.  B.  Beds,  »,  Rat. 

101  Uaahlnsts  Straat 

■awport,  Mwada  laland  1 

M-.  Daniel  BUabarg 
lha  RAID  Coryeratls 
1700  Main  Strst 

Santa  Monica,  Callfoada  1 

A1  aln  Bnthevs 
3009  Belly  Strst 

Alavandrla,  Tlrglala  1 

Profaaaer  Bbarhard  Ihla 
kk5  CL 

unlvaraity  ef  Plttabs^ 

Plttabs^  13,  ^nnaylvanla  1 


Bsas  of  Ssplluu  nod  Aceanta  Coda  OW 
Dapartsat  ef  tha  fcvy 
WMhl^tea  25,  D.  C. 


1 


Contract  Bsr*223(30) 


IteiTill  M.  Plood 

rnmtml  iMith  XMUtvt* 

20$  SorU  rer—t  4pih 
Imm  Artar, 


rnmtttf  0«l«l«r 
Ite  MD  Cer^eratleB 
1700  ItelB  0tr«*t 
Sant*  Mcjlea,  California 

FroTaanor  1.  Oalkwai 
Ilipartawt  of  Itethaantleo 
Vbl-wmnitj  of  MtaMota 
IMf  lyalln  Ih,  HlMMota 

Olanalal  Fomdanion  OrOMai^Wntor  Atadlnn 
OlaBBlal  rennAatlen  of  A^l.  fccaoaloo 
Italnaroltj  of  California 
Borkalajr  California 

M-.  K.  J.  D.  aiUlot 

Nnrtetlnf  Ranoaroh  Oopnrtaant 
Oils  Itetbitaoe  Cbanleal  Carp. 

27$  tflneboatar  Amna 
■•V  Bavnn  4,  Coonaetlott 

Dr.  M.  B.  Oomary 

IBI  Bn  ao  arch 

Torktown  Bal^ta,  York 

Profaaaor  2vl  Qrlllebta 
Departaeijt  of  ■eoocalet 
Onlvaraltj  of  Chleafo 
CbieaAO  yj,  ZUleola 

Dr.  laadora  Sailor 
■ary  tianaciarnt  Offlec 
WaahlACteo  2$,  D.  C. 

Er.  IbaodOTo  K«  Harrit 
fha  RAID  Corporatlco 
1700  Mala  Btraat 
Santa  Monica*  California 

Dr.  John  C.  Karaaajrl 
DopartMoat  of  Beonoaica 
Wnaraa  Btata  ttnl-reraltjr 
Oatrolt*  NlehlMan 

Profaaaor  M.  R.  Baotaaoa 
Oapartjwnt  of  Hatbaaatloa 
Onlvaralty  of  California 
boa  Ancolaa  2k,  California 

Profaaaor  C«  BllAratli 
NlchlMnn  Btata  Oolvartlty 
Boat  l^aalns,  Klcblaan 

Mr.  C.  J.  Bitch 
Iho  RAKD  Corporatlco 
1700  Main  Straat 
Santa  Mealea*  California 

Dr.  Alan  J.  BoffMno 
m  Baaaareh  Oantar 
Bn  218 

Torktova  Btlchta*  lev  York 

Prcfaaaor  Auatla  Boggatt 
School  of  Bualneao  Malnlatratloo 
Oolvaralty  of  California 
Btrkalap  California 

Ik.  C.  C.  Bolt 

Qrad  School  of  Induatrlal  ASalniatratlon 
ramagit  Inatltata  of  Vachaelogy 
PlctabwBb  13,  PaoBaflTaala 

Or.  John  V.  Bocqpar 
Coalaa  Powidatloa  for  Haaaarcb 
la  Bceoaalca 
Tala  Statloa*  Boa  2129 
B«v  Haraa,  Corawctlcut 

Pr^aaaor  B.  HetaUlng 
Dapt.  Of  IMthaantleal  StatAatlca 
thtlvaralty  of  Mortb  Carolina 
Chapal  Bill*  Bortb  Carollaa 

Profaaaor  I.  S.  Boathakhar 
Ikpartaant  of  Mcoaoalea 
Barrartf  Dalvaraltjr 
CaabrldAB  36*  l^aaarhMaatta 

Dr.  H.  M.  Bi^wa 
Bapartaant  of  Ucaatrlea 
School  of  Arlatloa  IMdielaa 
0.BJL.P. 

n«iA,  TmjOB 


iiSiaint  Belaaooa  Maa.  Pr«J. 
Liaraltar  «f  OhUforala 
I  M^laa  2k,  California 


Oapartaaat  of  Beoatnloa 
Bhlaaral^  of  Oalifnvla 
BvMMliar  k,  California 

JcMm  InplilBO  BMlnaraltf  Uhraiy 
Arf  laltlaaa  Bapai  taaai 
•iltlaara  lA,  Ifcrylaait 

Copt.  Mhltar  1.  Eaaa,  «■ 

Offlca  of  Baaal  Raaaareh 
CoAa  kOO 

Bapartaant  of  Baay 
tMahlnctoa  23,  O.  C. 

Profaaaor  taarraaca  Dala 
Dapartaaat  of  BcOBOalea 
Onlaaralty  of  Panaaplaanla 
FhUadalphla,  Paaa^laonla 

Profaaaor  T.  C.  YniTpaanr 
Cotflaa  Poaadatloa  for  Baaaareh 
la  Bcoaoalca 
Tala  Station,  Bn  212$ 

Bav  flaran,  CoBaaetleiit 

Profaaaor  Harold  M.  gv>>*i 
Bapartaant  of  Mathaaatlea 
Irpn  Maar  Collapa 
Myn  IMar,  Panaarlvaala 

Ik.  Boward  laitln 
lha  SAID  Corporatloa 
1700  Main  Straat 
Snta  Monica,  California 

Sr.  R.  P«  laaaallottl 
Bapartaant  of  Beonoalea 
WaahlBgtoa  Stata  CoUapa 
Pttllaan,  Nhahlnctoe 

Profaaaor  C.  B.  Laaka 
Dapartaaat  of  Mathaaatlea 
fitoaaalaar  Polpiachnlcal  Xaatltuta 
fr«p,  Bav  York 

Frofaaacr  W.  V.  taontlaf 
Dapartaaat  of  Bconealea 
Harrard  Uhlaaralty 
Ci^brldpa  36,  Maaaaebaaaita 


Canter  for  Intamatlonal  ATfalm 
6  DlTlnlty  Aaa. 

Inrmrd  Qnlvaralty 
Caakrldpa  3S,  IMaaaclanatta 

Profaaaor  8.  B.  Zdttaaar 
Dapt.  of  Indvatrlal  AigiBaarli^ 
Col^^la  Ualwralty 
Bav  Ycrk  27#  Rev  York 

Dr.  I.  Oncan  lura 
Dept,  of  Papcholecf 
Unirvraltjr  of  ^onaylvnala 
Phlladalphla,  Pana^lmala 

Sr.  Cralp  A.  Tjagirtri 
299k  Vlnehaatar  khy 
landko  Cordova,  Ckilfnnla 

lanolt  IMadalhret 
Xatamatloaal  Baalaasa 
IMchiaao  Corp. 

Haaaarcb  Cantor 

P.  0.  Bn  21B,  btata 

Terktown  Baipbta,  Mav  York 

Prcfaaacr  Jallva  Marpolla 
Dapnrtaaat  of  Baalnana  Malalatratlea 
Dniwralty  of  Callfomla 
Berkalap  k,  Califorala 

Ik.  Barry  M.  Hvkowlt* 

Loplatlea  Bapartaant 
Yka  lABO  Corporatloa 
1780  Main  Ptroat 
siata  Monica,  California 

Profaaaor  Jaeoh  tMrarliak 
School  of  BBafaaaa  dSaialfftratlan 
Ohlvoralty  «f  Califorala 
X<M  Aapalaa  2k,  Califorala 


trheol  af  BMinnna  ASaialatraUa 
Bhivarai^  of  Califorala 
larhalay  k,  CMUforaia 

Pwfaaaor  LAeaal  R.  tMBaaala 
Bapartanat  of  Beoaoslca 
Bhivaraity  af  Boekaatar 
Uvar  ttatloa 

Boebatar  20,  Bmt  York 


laad,  Oparatloaa  Baaaareh  BaaR 

V.  S.  B.  balalap  Oavloa  Onatar 

Port  VMaRl^toa,  L.  I.,  Bav  Tark  1 

Ik.  Biebrd  A.  Millar 
kcrri  Mnat  7U  itraat 

Port  Worth  7»  bxaa  1 

Profkaaor  Prance  NedlsUaai 
Qrad.  School  of  Indvatrlal  Adolalatratloa 
Paraacla  laatltvta  of  bohaolosr 
Plttahvrsh  13,  Naa^lvaala  1 

Prefcoaer  Jaaoo  Neora 

Dipt,  of  Zadvatrlal  Baplniarias 

lorthoaatom  Ohlvoralty 

Beaton  12,  Mm*.  1 

Profoojor  0.  NorpoBStan 
Bcencootrle  Baooar^  Prepr^ 

Prlacatea  Bhlvarsity 
92*A  Baaaan  Straat 

Prlneoton,  Bav  Jaraay  1 

D.  I.  Banhao 

Chlaf,  Indvatrlal  tiplnaari:^  Mv. 
C^ptroUar 

Bdarn,  Boa  Bamardlao  Air  btarlal  Area 
Borten  Air  Poraa  baa,  Califorala  1 

Ik.  M.  L.  Berdan 
Baaanreh  Mviiien 
CoUapa  of  Miplnaarlns 
Bav  York  Uhlvaralty 

Bav  York  $3,  Bav  York  1 

Prefoaaer  B.  X.  O'BaUl 
Sapartaant  of  BaslBaarlas 
Uhlvaralty  of  California 

loa  Anpaloh  2k,  California  1 

Ik.  Richard  T.  Palamo 
Oparatloaa  Hcaaareh  Bapartawnt 
VUlev  Rva  laheratorlao 

Tpallantl,  Mlehlpua  1 

Profaaaor  Stanley  Boltar 
Bepartaant  of  Icoaokleo 
Pvdva  Oaiwvolty 

lafayatta,  laaiana  1 

Profaaaor  M.  BUay 

301  Ptaphaaa  Hrrlal  bll 

Vhlvaralty  of  California 

Barkalay  k,  Califorala  1 

Preftaacr  B.  Booanblatt 

Dapartaaat  of  Ptatlotleo 

Itw  Oaorpa  Uhahlaptoa  Ohlvomlty 

kbhinptcn  7,  D.  C.  1 

Profooocr  Mkarior  Boaoahlatt 
Bopntaaat  of  Ikthaaatleo 
Bmn  Utolvaroity 

Previdaaco  12,  IBicdo  ^ 

Prefaoaer  A.  B.  Boos,  laod 

Oaportoaat  of  Nnttaaaatlea 

Oalwroity  of  lotro  Dhb 

Botra  Dona,  1 

Profaoaor  Jaroac  Bathabarp 

Dapnrtapnt  of  Beennlee 

Berthoaatam  Ohlvarolty 

Bvanatoa,  ZUlaoio  1 

Profoooer  Albart  I.  Rubanatala 
Dapartaant  of  XaSaatrlal  ^plniTlng 
Borthmatam  Dnivaralty 

Bmaatoa,  XUlaola  1 

ProfOaaor  Be  man  BWbla 

Dapartaaat  of  Statlatlea 

Mtchipin  Stata  Ohlvoralty 

Baat  laaalas,  Mlehl^a  1 

Profaaaor  Rl^ord  S.  Baiaor 

Daparbaat  of  fkUooophp 

Mtchlpna  Stata  Shlwralty 

Boat  baatap,  Wcht^a  1 


Coatract  Boat -229(30) 


Praf«*«er  P.  A. 

THyrTiiit  at  BeoMBlet 
M.  I.  T. 

CMlri  flt  y9,  IteaiarlHMitU 

Ik*.  Z.  RleharA 
Bchocl  of  BuImm 
▼lacMt  IkU 
VblToraitjr  of  Mlanaaot* 

HlBaMOta 

ProfoMor  4ndr«v  Setaati,  Jr> 

Dapt*  of  Indwtrlal  Incliwtfinc 
Coniell  Oalvoralty 
Ithaca,  lav  Tork 

Profaaaor  B.  A.  81aoe 
Dapt.  of  laduatrlal  Adatalatratieo 
Caraacl*  leatltttta  of  Taehnolocp 
PlttMwrgh  13»  PanaajrlTaala 

Ifr.  Jm  R.  8iiq>aon 

Buraau  of  Bi^liaa  aad  i^eomta 

Coda  as 

Arllactoe  Aonex 

U.  8.  Sapartaant  of  Ravy 

Waahlncton  25,  G. 

Profeaaar  laurle  Snell 
Departaest  of  NatheMtlca 
Dartmouth  College 
Banover,  Rev  Basgiohire 

Profeaaar  Robert  Solow 
Dapartmeot  of  BcoBCMlea 
Maaa.  loat.  of  Vaehnologjr 
Canbrldge  *0,  Maacachaaetta 

Profeaaor  Robert  Simaera 
Department  of  IcocKKlca 
Ublveraltj  of  Peonasrlvania 
Rilladelphla,  Pennaylvanla 

Profeaaor  Oerald  ‘Ataqpaan 
Dapt.  of  laduatrlal  Alalniatratioe 
Carnegie  Inatltute  of  Technology 
Plttabur^  13,  Pennajrlvania 

Profeaaor  R.  M.  Thrall 
Department  of  Matbamatlca 
tlnlvaraltjr  of  Mlehl«aB 
Abo  Arbor,  Xlchl^e 

Profeaaar  L.  M.  Tlehrinaky 
DaparteMnt  of  Rnglneerlng 
INilveral^  of  California 
Berkeley  California 

Profeaaor  Jamaa  Tobin 
Council  of  Iconoilc  Adviaers 
Rxecutlre  Office  of  The  Prtaldent 
Vaahlagton  2$,  D.  C. 


Profeaaor  C.  B.  Ttm^klna 
Dept,  of  Mathomatlce 
UolTvralty  of  Callfomla 
Loo  Angelee  2>»,  Callf<ml> 

Profeaaor  A.  V.  Tucker 
Pine  Hall,  Box  7C6 
Princeton,  Rev  Jeraey 

Profeaaor  Rlrofvd  Ueava 
Behavlwal  Seleacea 
202  J'nnlpero  Serm 
Stmnford,  California 

Profeeaor  0.  F.  Totmv,  Jr. 
Department  of  Icoeomiea 
Tale  Uulreralty 
Rev  Haven,  Cacaectlcut 

Profeaaar  V.  A.  UalUa 
207  HaakaU  BaU 
Ohlveralty  of  Chlcoge 
Chicago  37,  lUlnola 

Profeaaor  J.  L.  Waleh 
Department  of  Nathamvlca 
Harvard  OnlTaratty 
CMbrldga  Maeoachuoatta 

Or.  T.  Ubltin 

School  of  laduatrlal  MaoagMent 
Maaa.  Inatltute  of  Tk-chnolagy 
Ca^trld^  Maasaebueetta 

M*.  Meyw  A.  Vlckeljpw 
6393  •  A  Oaim  Stroot 
Oakland  9,  California 


Ife-.  ftalUp  Uolft 
Bw  BMB)  Cofparatioa 
1700  Mils  Street 

Baata  Monlea,  Ckllfcroia  1 

Pi'ofeaaor  J*  Velfnrita 
Tlufai  tmant  ot  MntlMHatles 
Ceraall  Vkaluoreltp 

Xtlimea,  Itow  Tet^  I 

Profeaaor  Robert  J.  Walfeoa 
CES  Ix»e. 

9171  Wilahlre  Blvd. 

Rerarljr  HlUb,  California  1 

Hr*  Iferohall  X.  Wood,  Chief 
Zi^natrlmX  fulaamhllltp  hrmaeh 
Offlea  of  Aaa't.  See.  of  Scfanae 
VaohlAgteax  2$,  D.  C.  1 

Profeaaor  M.  A.  Voodbwy 
Departamnt  of  Itethomatlea 
Rev  Tork  Ihilvnraltx 

Rev  Tork  53,  >•»  York  1 

Diatrlbutlon  via  Ot  London 

CcooMndlng  Officer 
B-aneh  Office 
Rory  Ro.  100 
Fleet  Poet  Office 
Rev  Tork,  Rev  Tcrk 

Profeeeor  Maurice  AUale 

62  BottloTard  St.  Nlehel 

Parla  6a,  FIURCE  1 

Mr.  William  N.  O^md 

UDlreralty  of  Blxmln^um 

RLmlngham,  *""Tf^»rn  X 

Mr.  Frank  Hahn 
Department  of  Icomcmlea 
nntrchlU.  Cellega 

CMbrld«e,  BRUARD  1 

Dr.  Rudotlf  Buaoer 

Inatltut  aagawBBdte  Ikthomatlk 

Dnlraralt&t  Bam 

Bam,  SWimiARD  1 

M.  Jaeqoea  laootme 
SoeldtA  dme  Ratbamatlguas 
Appll«wfoa 
33  rue  Matador 

Parle  9»»  FMRCR  1 

Dr.  Rrtmoewl  Mallnvaud 
Inatltute ^tetleoal  de  la  Stattatlque 
at  dea  Rtudea  Xccmcailguae 
?9,  Rual  Bnaoly 

Parle,  FRARCX  1 

^ofeesor  Meurice  McNanue 

Facttltgr  of  Coemeree 

and  Social  SeloDce 

Itelwrelty  of  BlmlnghM 

Blruln^iaai,  Rngland  1 

Frofesaor  Dr.  Balcc  Bauamaan 
Inatltut  fvr  Prandeovurkehrevlaaenaehaft 
Orefatrmeae  39 

Frankfurt  a.M.,  ODBMRT  1 

Profeaeor  Henri  Ihell 
Bcooometrlc  Inetltute 
■etberlmnde  Bdioel  of  Rccocmlea 
Rotterdem,  RRPHBUARDB  1 


V.  6.  Attache' 

Attache'  for  Air  (Tokyo,  Japan) 
APO  •  5CO 

Ban  Praocleco,  California 

Frofeeeer  Nlcblo  Norlahlma 
Inatltute  of  Boelal  and  Rceocmle 
Re  near  eh 

Ublvuralty 
Tuyooaka,  Oeaka 
JAPA 

Tekaahl  Reglahl 

112  Barvmeho  Magvoku 
Tokyo,  JAPW 

Dr.  —  Rikalde 

Zaatltute  of  Boelal  and  Bcomtale 

Reeaarcb 

Ihilvurolty 


Refartaest  ot  leeM 
Ttatefea  Bklvurslty 
Saadai,  JAPAR 


Jan  Ikenta 

of  BeaBaml.ea 

Dtaiueraity  of  R.  B.  Uklea 
Kanalngtoa,  R.  8.  Nhlaa,  ADBOULZA 

FTofoaeor  David  Slater 
Qeeen'a  Rtalvoralty 
Kingetom,  Oetario 


Profeeaor  D.  T.  T.  Baer 
Dapt.  Seonomlca 
Vnlvereity  of  Chicago 
Chicago  37,  ruinola 

Addltlcnal  ccplaa  for  Rrejeet  leader 
and  aaelataBts  and  rasarvu  for  future 
ra^ulromenta . 


Cootraet  Rear»223(30) 
BiiTveaber  196I 

(2a) 


